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£\j , We present a perturbative calculation of the neutron matter equation of state based 

on low-momentum two- and three-nucleon interactions. Our results are compared to the 
model-independent virial equation of state and to variational calculations, and we provide 
theoretical error estimates by varying the cutoff used to regulate nuclear interactions. In ad- 
dition, we study the dependence of the BCS 1 So superfluid pairing gap on nuclear interactions 
and on the cutoff. The resulting gaps are well constrained by the nucleon-nucleon scattering 
phase shifts, and the cutoff dependence is very weak for sharp or sufficiently narrow smooth 
regulators with cutoffs A > 1.6 fm -1 . 

§1. Introduction 

o i 

The determination of a reliable equation of state of nucleonic matter plays a 
central role for the physics of neutron stars 1 ) and core-collapse supernovae. 2 )' 3 ) Fur- 
thermore the superfluidity and superconductivity of neutrons and protons is an im- 
portant phenomenon in nuclear many-body systems, 4 - 1 ' 5 ) in particular for the cooling 
t^j- ' of neutron stars. 6 ) In this contribution, we present calculations of the neutron mat- 

ter equation of state at finite temperature and of the 1 So superfluid gap in the BCS 
approximation based on low-momentum interactions. 

Renormalization group methods coupled with effective field theory (EFT) offer 
the possibility for a systematic approach to the equation of state. By evolving 
I nuclear forces to low-momentum interactions Mowfc 7 ^ 9 ^ with cutoffs around 2fm~ 1 , 

the model-dependent short-range repulsion is integrated out and the resulting low- 
momentum interactions are well constrained by the nucleon-nucleon (NN) scattering 
data. Furthermore, the corresponding leading-order three-nucleon (3N) interactions 
(based on chiral EFT) become perturbative in light nuclei for A < 2fm~ 1 . 10 ) 

With increasing density, Pauli blocking eliminates the shallow two-nucleon bound 
and nearly-bound states, and the contribution of the particle-particle channel to bulk 
properties becomes perturbative in nuclear matter. 8 ) The Hartree-Fock approxima- 
tion is then a good starting point for many-body calculations with low-momentum 
NN and 3N interactions, and perturbation theory (in the sense of a loop expansion) 
around the Hartree-Fock energy converges at moderate densities. This can be under- 
stood quantitatively based on the behavior of the Weinberg eigenvalues as a function 
of the cutoff and density. 8 )' 9 ) 
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Some uncertainty remained concerning a possible dependence of the So pairing 
gap on the input NN interaction in low-density neutron matter (hp < 1.6 fm ). 
We address this point and explore the dependence of 1 So superfluidity on nuclear 
interactions at the BCS level in detail. We find that the BCS gap is well constrained 
by the NN phase shifts. Therefore, any uncertainties are due to polarization (induced 
interaction), dispersion and three- nucleon interaction effects. 

§2. Equation of State of Neutron Matter 

Using the Kohn-Luttinger-Ward theorem, n )> 12 ) the perturbative expansion of 
the free energy (at finite temperature) can be formulated as a loop expansion around 
the Hartree-Fock (HF) energy. In this work, we include the first-order NN and 3N 
contributions, as well as normal and anomalous second-order NN diagrams. Other 
thermodynamic quantities are computed using standard thermodynamic relations. 
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Fig. 1. Energy per particle E/N as a function of the density p at first order (left panel) and 
including second-order NN contributions (right panel). 13 ' 

The resulting energy per particle E/N as a function of the density p is shown 
in Fig. 1 for a cutoff A = 2.1 fm -1 and temperatures T = 3, 6 and lOMeV. 13 ^ The 
results presented in the left panel are the first-order NN and 3N contributions, and 
those in the right panel includes all second-order diagrams with NN interactions. 
For T = 6MeV, we also give a band spanned between A = 1.9 fm" 1 (lower line) and 
A = 2.5 fm -1 (upper line). The inclusion of second-order contributions significantly 
reduces the cutoff dependence of the results. The model- independent virial equation 
of state 14 ) and the variational calculations of Friedman and Pandharipande (FP) 15 ) 
are displayed for comparison. 

The inclusion of second-order correlations lowers the energy below the variational 
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Fig. 2. The neutron-neutron 1 So superfluid pairing gap A as a function of the cutoff A for three 
densities and different smooth exponential regulators, as well as for a sharp cutoff. 18 ' The 
low- momentum interactions are derived from the N 3 LO chiral potential of Ref. 19 ' 



results for densities p < 0.05 fm~ 3 , and we observe a good agreement for E/N with 
the T = lOMeV virial result when the second-order contributions are included. In 
the virial equation of state these contributions are included via the second-order 
virial coefficient, while in the variational calculation the state dependence of such 
correlations is only partly accounted for. 16 ) Furthermore, the generic enhancement 
of the effective mass at the Fermi surface leads to an enhancement of the entropy at 
low temperatures above the variational and HF results. 13 ^' 16 )' 17 ) 



§3. BCS gap in the 1 Sq channel 



We solve the BCS gap equation in the 1 So channel 

A(k) = -i / dpp* Y^MM , (3-1) 

with the (free-space) low-momentum NN interaction Vi ow k(k, k'). Here = e(p) — 
fi, e{p) = p 2 /2 and fi = k\j2 (c = h = m = 1). 

We find that the neutron-neutron BCS gap is practically independent of the NN 
interaction. 18 ^ Consequently, 1 So superfluidity is strongly constrained by the NN 
scattering phase shifts. The maximal gap at the BCS level is A « 2.9 — 3.0 MeV for 
kp ~ 0.8 — 0.9 fm -1 . For the neutron-proton 1 So case, we find somewhat larger gaps, 
reflecting the charge dependence of realistic nuclear interactions. 18 ) 

In Fig. 2 we show the dependence of the neutron-neutron 1 So superfluid pair- 
ing gap on the cutoff starting from the N 3 LO chiral potential of Ref. 19 ) for three 
representative densities. 18 ^ We employed different smooth exponential regulators 
f(k) = exp[— (k 2 /A 2 ) n ], as well as a sharp cutoff. As long as the cutoff is large com- 
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pared to the dominant momentum components of the bound state (A > 1.2/cp), the 
gap depends very weakly on the cutoff. This shows that the 1 So superfluid pairing 
gap probes low-momentum physics. Below this scale, which depends on the density 
and the smoothness of the regulator, the gap decreases, since the relevant momentum 
components of the Cooper pair are then partly integrated out. 

§4. Conclusions 

In summary, we have studied the equation of state at finite temperature including 
many-body contributions in a systematic approach. We have found good agreement 
with the virial equation of state in the low-density-high-temperature regime. Ana- 
lyzing the cutoff dependence of our results provides lower bounds for the theoretical 
uncertainties. The possibility of estimating theoretical errors plays an important 
role for reliable extrapolations to the extreme conditions reached in astrophysics. 

In addition, we have shown that the 1 So superfluid pairing gap in the BCS 
approximation is practically independent of the choice of NN interaction, and there- 
fore well constrained by the NN scattering data. This includes a very weak cutoff 
dependence with low-momentum interactions V\ ow k for sharp or sufficiently narrow 
smooth regulators with A > 1.6 fm -1 . At lower densities, it is possible to lower the 
cutoff further to A > 1.2kp. Furthermore, the pairing gap clearly reflects the charge 
dependence of nuclear interactions. The weak cutoff dependence indicates that, in 
the 1 So channel, the contribution of 3N interactions is small at the BCS level. 
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